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We design composite controlled-phase gates, which compensate errors in the phase of a single gate. The 
errors can be of various nature, such as relative, absolute or both. We present composite sequences which are 
robust to relative errors up to the 6th order, with the number of the constituent gates growing just linearly with 
the desired accuracy, and we describe a method to achieve even higher accuracy. We show that the absolute 
error can be canceled entirely with only two gates. We describe an ion-trap implementation of our composite 
gates, in which simultaneous cancellation of the error in both the pulse area and the detuning is achieved. 
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I. INTRODUCTION 

Two-qubit gates form the basis of quantum computing: any 
quantum computation can be constructed entirely by these 
Mtes, combined with the one-qubit Hadamard and phase gates 
lUn . Prominent examples are the controlled-phase (CPHASE) 
gate and the closely related controlled-NOT (CNOT) gate, 
which have been implemented on various physical platforms, 
such as ion traps lH , nuclear magnetic resonance 1^, photonic 
qubits 101, superconducting qubits Ist] and atoms . 

Quantum computing, however, depends critically on the 
accuracy with which these gates are implemented and fault- 
tolerant computation is possible only if gate infidelity is very 
low, typically below 10 §]. Different sources of error 

can be identified in an experiment ranging from decoher¬ 
ence to imperfections in the control fields, most notably pulse 
length errors, field inhomogeneity, improper gate duration, 
frequency shifts, etc. Because errors in the control fields lead 
to incorrect rotation angle of the state vector on the Bloch 
sphere, we refer to them as rotation errors. 

Various techniques to deal with rotation errors have been 
proposed and implemented. Composite pulses, for example, 
have found broad application, where systematic errors in the 
control field play a major role (1110]. A composite pulse is 
a sequence of pulses with well-defined phases, which are de¬ 
signed such that the errors from the constituent pulses largely 
cancel each other. For example, the very popular broadband 
composite pulse BB1 of Wimperis |0] cancels the error up to 
the second order on arbitrary state (without assumption of the 
initial state of the system). More accurate pulses have been 
derived by Brown et al. lU^ but they have found limited ap¬ 
plication as they quickly become extremely long. Low et al. 
Illin found optimal sequences yielding the same accuracy as 
Browns’ but with much fewer primitive pulses. 

In an important development, Jones ifl^ showed how the 
available single-qubit composite pulses can be used to con¬ 
struct composite conditional two-qubit gates, which are the 
backbone of quantum computation. In particular, he extended 
the BB 1 pulse by Wimperis to construct a second-order broad¬ 
band two-qubit gate, which is robust to rotation errors. Later, 
Hill foil showed how to achieve robust CNOT gates from al¬ 
most any interaction based on BBl and used sequence con¬ 
catenation for higher precision. Gates of higher precision 
can be obtained even more efficiently by extending Low’s se¬ 


quences d. 

In this work, we derive improved highly-accurate compos¬ 
ite CPHASE gates of shorter duration and length than pro¬ 
posed so far. Although derived for two qubits, these sequences 
can be used to construct robust multiqubit CPHASE gates as 
well. In Sec.im we present sequences, which cancel relative 
errors up to the 6th order using up to 12 gates, and describe 
how to obtain even higher accuracy. In Sec. [Till we design 
sequences, which can handle errors of both relative and abso¬ 
lute nature. Then, in Sec. lIVl we describe an implementation 
of our composite gates with linear laser-driven ion traps. Re¬ 
markably, our gates can compensate simultaneous errors in 
the Rabi frequency and in the detuning. 


II. COMPENSATION OF RELATIVE ERROR 
A. General framework 

An ideal two-qubit CPHASE gate, denoted by (0), is rep¬ 
resented by the propagator (in a rotated basis) 

t/(0) (1) 

where Gx is the Pauli’s x matrix and 0 is a rotation angle Q. 
When two qubits are coupled with a coupling constant J for 
a time T, then 0 = JT. Relative rotation errors are described 
by multiplying 7 or T by an unknown factor, 1 -|- e, so that 
0 is higher or lower than a desired value 0. Therefore, rather 
than (0), in reality one obtains (0(1 -b e)), represented by the 
propagator t/(0(l -be)). The sensitivity of the gate to e can 
be much reduced by replacing single rotations with composed 
rotations, 

[<Pa'+i] {On) [^a] • • • (0i) [^’i] (®o) [^] ■ (2) 

Here time is running from right to left, so that the rightmost 
gate is the first one applied. This sequence contains V -b 1 
CPHASE gates (0^) and V-b2 single-qubit phase gates [tpi^], 
applied to a preselected qubit; once chosen, the same qubit is 
used over the entire sequence. We will apply [tpi^] on qubit 2. 
The sequence (|2]i is represented by the following propagator 

u^^\e) = F{(pN+i)uieN)F{(pN)---uieo)Fi(po), o) 
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where 

F{(p) = e-‘^^-~. (4) 

Note that by using the property (lAll) one can incorporate 
the phases in the U -gates, yielding 

[/W(0) (0l)t/^(0o), (5) 

or t/W(0) = ^’(^Af+i)nLo% in a more compact 
form. Here we have defined a phased CPHASE gate (0)^, 
represented by the propagator 

t/^(0)=e'®^■'^^ (6) 

where (J,|, = CJ^cos^ + qvsin^, Uo{9) = U{9) and (0)o = 
(9). The realization (|5]l may be more convenient in a practical 
setting, where the operator t/^ ( 0 ) is achieved at no additional 
cost by a simple shift of the phase of the driving field, which 
does not represent a physical modification of the qubit. We 
will use the realization Q for our composite sequences. The 
following relations can be derived from the property (lAlk 

l-l N 

0/ = -2 ^ (Pk, ^N+l = 52 (7) 

k=0 k=0 

with I = 1,2,... ,N. 

Three families of composite sequences are generally con¬ 
sidered; broadband, passband and narrowband, the first two 
of which will be of interest in the present work. 

B. Broadband sequences 


sequence (|5]l. It follows from Eq. (IA2b that {9k) is equal 
to 1 for 9k = n (even s) and to iaxOiji^ for 9k = j (odd s). Thus 
it can be shown that the product n^=i {9k) applies cj"’ to 
qubit one, where m counts the gates with 0^. = in the prod¬ 
uct. Because the zeroth-order approximation of B„ (0) must 
reproduce U (0), we must have H^i {9k) = 1 for e = 0. 
This implies that m is even, so that besides the target propa¬ 
gator U (0), the sequence Q must contain an even number of 
^-gates. Hence, the gate sequence (|5]l acquires the form 

W(f)^Jf)^,(©)o, forN = 2, (9a) 

forN>2, (9b) 

having length of A -I- 1 and a total angle of -f 0. 

We further found that in our sequences with A > 6 and 0 = 
j we can set = 0. The important implication from here is 
that we can merge the first two gates into a single gate (^). 
Note that, up to a global phase of n, this gate is equivalent 
to (1)^. As a result, we can eliminate one gate (^) from 
the sequences (l9bl i simply by setting 0o = thereby yielding 
shorter sequences, 

fo>-A>6, (10) 

with length A and a total angle of A^ — ^. We have subtracted 
1 from each index for consistence of notation. 

With the above assumptions the left-hand side of Eq. ([8]) 
can be handled relatively easy by applying the identities (IA31 l. 
(IBlb and (IB2b . shown in Appendices lAl and iBl Eor / = 0, we 
obtain 


I. General principles 


While every gate in the sequence (|5]l is first-order sensitive 
to e (all angles 9k are systematically wrong by some constant 
fraction e), the phases (j)k can be chosen such that the compos¬ 
ite gate sequence is robust to e up to a certain higher order n. 
To this end, we nullify the n lowest-order propagator deriva¬ 
tives with respect to e by solving the following system of n -f 1 
algebraic equations for the phases 0 *.: 


de‘ 


{/{A') (©)-{/(©) 


= 0 , 

£=0 


(8) 


^(— 1 )^ ^(j)k + (pN+l =0. 


k=0 


( 11 ) 


Higher-order terms are not simple enough to be useful. 


2 . Two-pulse sequence, n = 1 

Let us consider a sequence with A = 2. Eollowing the above 
arguments, we set 0i = 02 = f. Zero- and first-order errors 
are cancelled by imposing the following set of equations: 


with / = 0,1,...,«. Such sequences exhibit robust profiles 
vs £ around 0 and are called broadband. We denote them as 
B„(©) below. Longer sequences provide more free parame¬ 
ters to vary {9k and (j)k), thereby allowing to eliminate higher 
orders of £. 

Numerical calculations indicate that we must have 0o = © 
and 00 = 0 , and also that = 0 for all sequences with 
A > 2. Thus we are left with 2A — 2 parameters 9k and 0^ to 
solve for. It can be shown that for k>Q the angles 9k can take 

values n{s 1 /2), where s = 0,1,2,_ In what follows, we 

will restrict ourselves to 9k = ^ or n, in order to minimize the 
total angle and thereby the total time duration of the composite 


-01 + 02 + 03 = 0 , ( 12 a) 

f = 0 . ( 12 b) 

We obtain 0i = 0 , 02 = 30 and 03 = —20, where 0 = 
arccos(—©/;r), and the sequence is 

B,(©) = [-20](f)3^(f)^(©)o. (13) 

We consider this sequence to be of significant interest from ex¬ 
perimental viewpoint for its reasonable robustness and small 
duration and length. 
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B„(0) 

total angle 

phases 0o,01 , 02,03, • • ■, 

Bi(0) 

1.25;r 

0, 0, 30, —20 [with 0 = arccos(—0/;r)] 

B2(0) 

2.2571 

0, 0, 30, 0 [with 0 = arccos(—6/2;r)] 

B3(tt/4) 

3.25k 

0, 1.725, 0.244, 1.127, 0.351, 1.785, 1.042 

B4(tr/4) 

3.75k 

1, 0.170, 0.170, 1.374, 0.677, 1.598, 1.818, 

0.528, 1.995 

B5(^/4) 

A.75k 

1, 0.065, 2.257, 1.826, 1.020, 0.487, 1.452, 

1.671, 0.132, 0.812 

B6(jr/4) 

5.75k 

1, 2.193, 1.933, 0.737, 1.932, 1.286, 0.641, 

1.531, 1.983, 1.240, 2.077, 0.579 


Table I: Broadband composite sequences B„(0), which cancel the 
relative error e up to order n, cf. Eqs. ij. The phases for n > 2 are 
given in units of n. The fidelities are shown in Fig.[T] 


3. Four-pulse sequence, n = 2 

The composite gate with = 4 corrects for e up to the sec¬ 
ond order. The following sequence is obtained 

B2(©) = (f)^W3^(f)^(©)o (14) 

with 0 = arccos(—0/2;r). The length is reduced to four gates 
as two adjacent phases are found to be equal. This semence 
coincides with the BBl pulse derived by Wimperis and 
later used by Jones in to construct a robust two-qubit gate, 
as discussed in the Introduction. 


4. Higher sequences, n>3 




Relative error £ 


Figure 1: Fidelity of our broadband composite gates B„(0) to re¬ 
produce the gate exp (if CTiCTj), forming the basis of the CNOT gate, 
versus the relative error e. The order n is displayed on each curve; 
0 corresponds to a single uncorrected gate (cf. Fq. (d)). The dashed 
line represents the lO^'* benchmark level. Note the dramatic increase 
of the ultrahigh fidelity range (infidelity below 10^^) with N. 


Sequences of higher accuracy are calculated numerically 
(see AppendixICli. We picked a target angle of 0 = ;r/4, which 
is traditionally used in quantum information to construct the 
CNOT gate 111]. Third, fourth, fifth and sixth orders in e 
are eliminated for N = 6, 7, 9 and 11, respectively (cf. Eqs. 
(l9bl i and (fTOl i). with corresponding total angles of A = 3.25n, 
3.15n, 4.75n and 5.75n. The phases of the composite gates 
are given in Table|I] where in places adjacent gates have iden¬ 
tical phases. Therefore, as for B2(0), we can combine these 
gates and reduce the overall length of the sequences. To our 
knowledge, apart from B2(0), which coincides with BBl by 
Wimperis |^, all broadband gates are original. 

The corresponding fidelities F vs the error £ are shown in 
Figure [U where the standard definition of F is used. 


Tr(Al'B) 

Tr(AtA) 


^Tr(A^B) 


(15) 


with A = U (0) and B = B„(0). Note that, as usual, a prop¬ 
agator infidelity of order £^" corresponds to an error term in 
the underlying propagator of order £". A comparison with the 
benchmark of 1 — 10^^ (horizontal dashed line) reveals that a 
single gate can be fault-tolerant only if the error |£| does not 
exceed 1.8%. The N = 2 composite gate Bi(;r/4) exhibits a 
clear improvement over the uncorrected single gate: the tol¬ 
erance range is already |£| < 11%. As expected, the longer 


composite gates are fault-tolerant over a wider error range: 
we have |£| < 22% for B2(7r/4); |£| < 30% for B3(7r/4); 
|£| <37%forB4(7r/4); |£| < 42% forBsjTr/T); |£| <46%for 
B6(;r/4). Our gates B4(;r/4) and B6(;r/4) compare very well 
with the broadband gates BB4 and BB6 by Low et al. d: 
our B4(;r/4) with total angle 3.75n performs slightly better 
than Low’s BB4 with total angle 4.25;r and B6(;r/4) with to¬ 
tal angle 5.75n performs slightly better than Low’s BB6 with 
total angle 6.25n. 

The total angle, which determines the duration of the en¬ 
tire sequence, is iin + & (with 2n -f 1 gates) for n <3 and 
(« — 5 ) ;r -f- 0 (with 2n gates) for 3 < « < 7, to correct the prop¬ 
agator to order n. For example, B 3 (j), which is third-order 
insensitive, has a total angle of 37r + 0. Note for compari¬ 
son that the same performance is achieved by the pulse B4, 
derived by Brown et al. d and used by Jones ifl^ . which 
requires 29 single gates and a total angle of 407r -T 0. 


C. Passband sequences 

Our broadband sequences are useful for the implementa¬ 
tion of highly accurate gates on an isolated qubit pair. In an 
actual experiment, however, it is possible that neighbouring 
qubits are involved in the interaction too, against our will, e.g. 





















4 


P(!1,K2(0) 

phases tj>Q,(l)i,(j)2,th,---,<l>N 

Pi,i(0) 

0, (j), —(j) [with (p = arccos(—0/27r)] 

P2.l(0) 

0, -Xi . -Xi +Xi,X\+X2,X\-X2,-Xi-X2,'’t~Xi 

Pl.2(0) 

0, (Cl, X\ + 2 : 2 . -X\+X2, -X\ -X2, X\-X2,'’t + X\ 

P2.2(0) 

0, (j), —(j), —0, 0 [with 0 = arccos(—0/47r)] 

Pl.3(f) 

0, 0.076, 1.604, 1.851, 0.595, 1.443, 0.751, 0.691, 1.111 

P3,3(f) 

1, 0.091, 0.644, 1.866, 0.941, 1.596 


Table II: Passband sequences P„, „,(©), which cancel the error e 
up to order n\ around e = 0 and up to order 112 around e = — 1 
(cf. Eqs. l ll6b l. The phases for Pi 3(0) and P3,3(0) are given 
in units of 7t. The fidelities are shown in Fig. We have 

Xi = arccos + Xl = arccos (—^ 4^+02 )■ Note that 

P2,l (0) is obtained from Pq2(0) by a sign flip of Xi- 



- 1.0 -0.5 0.0 0.5 1.0 


Relative error s 

Figure 2: Fidelity of our passband composite gates Pni .n, (f) versus 
the relative error e. The orders (ni ,n2) are displayed on each curve; 
0 corresponds to a single uncorrected gate (cf. Fq. Q). 


as a result of residual laser light addressing these qubits. To 

suppress this effect while maintaining the robustness of the and for P3,3(f) (with total angle 5.75n) we have 
broadband sequences, one can use passband sequences. 

Passband pulses satisfy the equations 


d'i 

de'i 

de‘2 


{/(A') (@) _[/(©) 

[/W(0)-ll 


= 0 , 

£=0 



(16a) 

(16b) 


where /i = 0,1,...,ui and I 2 = 0,1,...,n 2 - Equations ( I16ab 
define the broadband part around e = 0, and Eqs. (I16bl) de¬ 
fine the narrowband part at e = — 1; the latter ensure that 
small coupling strengths, as “felt” by neighbouring qubits, 
yield negligible rotation. These passband pulses, denoted as 
P„i „2(0), are robust up to order ni around £ = 0, and up to 
order U 2 around £ = — 1. As such, passband sequences real¬ 
ize robust rotations upon our pair of qubits, as achieved using 
broadband sequences, while suppressing rotations upon the 
remaining qubits S. 

Eqs. (fT^ can be handled relatively easy by applying the 
identities (IA3b . (IBlb and (IB2b . Eor Eqs. (I16ab we proceed 
as for the broadband sequences. For I 2 = 0, Eqs. (I16bb are 
automatically fulfilled, while for (2 = 1 and I 2 = 2, they are 
reduced to 


N 

20-t-7r^e'‘^<^ =0, (17a) 

A = 1 
N 

37t^-2&^ + 7t^ ^ = 0, (17b) 


respectively, which are treated numerically, as it is the case 
also for I 2 > 2. 

Calculated passband sequences correcting to different or¬ 
ders III and U 2 are shown in Table. El The sequences are as 
follows; for Pi 2 ( 0 ), P 2 .i(©) and Pi. 3 (f) (with total angles 
3n + &, 371-1-0 and 4.257r) we use the sequence ( l9bb . for 
Pi. 1 ( 0 ) and P 2 , 2 (©) (with total angles 27r-|-0 and 47r-l-0) 
we have 


(18) 


The corresponding fidelities for 0 = ^ are shown in Fig. |2] 
Some of our sequences can be found in Ref. ill: Pi,i(©) 
is identical to API, P2,2(0) is identical to PD2 and PBl by 
Wimperis Jgt]. Note that P 3 , 3 (f) with total angle 5.757r per¬ 
forms almost as good as AP3 and PD4 with respective total 
angles 6.257r and 8.257r. 


III. COMPENSATION OF ABSOLUTE ERROR 

In addition to the above sequences, which compensate rel¬ 
ative errors in the target rotation angle 0, we have designed 
composite sequences, which suppress absolute errors that oc¬ 
cur as a constant offset ^ in the rotation angles, {&)^ 

(0 + ^)^, represented by (0-f ^). Like relative errors £, 
these errors must enter systematically in the sequences. 

Using the property (/;[ (0) = Uq (—0), we have found that 
absolute errors can be eliminated completely from (0) with 
the sequence 

(®)A.^ = (-f).+4f)^^ (20) 

represented by 

t/A.0(0)=t/;i+0(-f)t/0(f). (21) 

Indeed, we have 

Ua,^ (0 + ^) = t/,+0(-| + ^)t/4| + ^) = 

+ ^)=t/A,^(0). (22) 

Potential errors in the phase ^an be removed following the 
composite technique of Ref. 11^ . 

Finally, we construct composite gates robust to errors of ei¬ 
ther nature, relative and absolute. This is done by substituting 
(0) with Ua,(I> (0) throughout in our sequences in Sec. HJ 


(®)o 
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For example, a gate robust to e to the third order and to ^ to 
any order is obtained with 

[-20](?)a,30(i)a.0(®)a,o> (23) 

where (j) = arccos(—©/;r). 

Below we discuss a physical realization of our composite 
CPHASE gate with laser-driven linear ion traps. While we 
consider ion traps, we note that our sequences are applicable 
to other systems, as well. 


IV. IMPLEMENTATION WITH TRAPPED IONS 

In trapped ions a pop ular two-qubit gate is the Sprensen- 
Mplmer (SM) gate IIO . It has been used by numerous ion 
trapping groups as a paradigmatic gate for quantum informa¬ 
tion processing. The SM gate was demonstrated by Leibfried 
et al. ifT^ and later by Kirchmair et al. ifl^ with fideli¬ 
ties around 97%. A record gate fidelity of 99.3% has been 
achieved by Benhelm et al. 11201] . Various dynamical decou¬ 
pling techniques have been experimentally demonstrated to 
protect the two-qubit gate from the environment j^ . Coher¬ 
ent error suppression using a pulse shaping technique has also 
been experimentally demonstrated 11^ . where the effects of 
certain frequency and timing errors were suppressed. 

In this section we will show how one can achieve a compos¬ 
ite CPHASE gate, robust to i) rotation errors, which may be 
caused by improper laser intensity and timing, and ii) certain 
frequency errors, which may be caused by a shift in the trap¬ 
ping frequency. The gate duration grows only linearly with the 
achieved precision contrary to previous proposals for coherent 
error suppression ll^ . where exponential scaling is observed. 


A. Hamiltonian and propagator 

Consider two ions irradiated along the transverse x direc¬ 
tion with a bichromatic laser field with frequencies otr = 
~ ®cm — ^ and tUb = Wo + Wem + tuned close to the 
first red (Wr) and blue (Wb) sidebands of a common vibrational 
mode. Here Wq is the frequency of the internal atomic transi¬ 
tion of each ion, Wcm is the frequency of the vibrational mode, 
and A is a suitably chosen detuning. The laser frequencies Wr 
and Wb sum up to twice the qubit transition frequency, while 
neither of the lasers is resonant to any level. Thereby only 
transitions where the atomic states are changed collectively 
take place. The interaction Hamiltonian is 

^ ^ E ) , (24) 

<:=1 ^ ^ 

where g is the (time-independent) Rabi frequency of the spin- 
phonon coupling and <t(C^^) = -p , with <7^ 

(CT^) being the spin raising (lowering) operator for ion k. The 
spin and the motional laser phases are defined hy = 

Q), where and are, respectively, the laser phases of the 
blue- and red-detuned laser beams as seen by ion k. 


The propagator U is obtained using the Magnus expansion 

iH: 


U =D(a)exp 



(Ar-sinA7')(j(Ci+)(T(C2+) 


(25) 


where T denotes the duration of interaction. D{a) is a dis¬ 
placement operator, D{a) = exp (aa^ — d^a) with 


a = - 


SJL 

A 




E 

k=\ 


y 




(26) 


which causes an undesired change of the vibrational state of 
the ion system. 

Now we discuss how to restore the vibrational state (elimi¬ 
nate D{a) from the propagator U), while preserving the con¬ 
ditional dynamics, described by the a-(J term in U. Note that 
if we only shift with n in Eq. (|26] |, we get a displacement 
of opposite magnitude, D{a) D{—a). This phase shift can 
be achieved either by a direct manipulation of the laser phase 
or by sandwiching U with 7t pulses on both ions. Therefore, 
in order to restore the vibrational state, we apply a second 
bichromatic pulse of equal Rabi frequency g and duration T 
with a phase shifted with n. Then the propagator becomes 


U = exp 


■.V 

' A2 


(A7'-sinAr)cj(Ci+)cT(C2+) 


(27) 


provided that potential errors in the interaction variables are 
systematic. Note that hereby we restore the vibrational state 
without even knowing the exact size of the detuning A. 

In the rest of the section we consider the implementation of 
both realizations Q and Q. The first allows us to use global 
addressing, a key advantage in the Sprensen-Mplmer gate 
id, at the expense of additional single-qubit phase gates, 
while the second requires individual addressing and possibly 
just a single phase gate (for N = 2). 


B. Implementation with global addressing 


Global addressing implies equal laser phases for both ions, 
i.e. = 1 ^*. Without loss of generality, we can assume that 
= 0, which implies that C7(i^^) = The propagator 
becomes 


U = exp 


4c ^ 

(^^ “ sinAT) 


(28) 


which yields the gate (cf. Eq. ([T])) with 0 given by 

4^2 

0 = -^ (AT-sin AT). (29) 

An important implication follows from here; potential sys¬ 
tematic errors in the Rabi frequency g (including unequal cou¬ 
plings), the detuning A and the pulse duration T combine into 
a single error in the rotation angle, 0 —> 0(1 -f e), which we 
already know how to suppress by using our composite broad¬ 
band sequences, listed in Table |T] 
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When there are more than two ions in the trap, residual laser 
light is likely to couple neighbour ions, as well; neighbour ion 
k will be coupled with Rabi frequency gk. As a result a ro¬ 
tation will occur with small angle 0*. = (Ar — sinAr), 
where we expect that 9^ <C 0. This effect can be well sup¬ 
pressed by using our passband sequences, listed in Table HIl 


C. Implementation with individual addressing 


Now the spin phase ^2 of ion two is modulated (relative 
to ion one), where the goal is to absorb the phase gates in the 
rotations [cf. sequences Q and (|5])]. Again, we set 
and without loss of generality, we assume that = 0. The 
propagator (l27l i becomes 


U = exp 


/^(A7'-sinAr)(J,,i(T(C2+) 


(30) 


The following identities are useful for calculating the error 
terms 


21 


=exp /£(-1)^0^(7; 


A:=l 

2/+1 


k=l 

21+1 




A:=l 


A=1 


(A3a) 

(A3b) 


where I = 0,1,2,.... For / = 1 we have 

(j(0i)( 7(02) =exp[-!(0i -(/) 2 )(tJ (A4a) 

and 

(7(0i)(7((/)2)c 7(03) = (T(01 -02 + 03). (A4b) 


which yields the gate [cf. Eq. ©I that we need 

for the sequence © with 0 given by Eq. ( |29] |. 


V. CONCLUSION 

We have derived highly-accurate broadband and passband 
CPHASE gates, which correct rotation angle errors of relative 
and absolute nature. Eor relative errors, the number of the in¬ 
gredient gates and the duration of our sequences grow linearly 
with the leading error order, as opposed to most proposals, 
where exponential growth is observed. Absolute errors can be 
eliminated completely with a sequence of just two gates. Im¬ 
plementation with trapped ions using bichromatic laser fields 
is discussed, where our sequences compensate errors both in 
the pulse area and the detuning. 


Appendix B: Calculation of propagator derivatives 


To calculate the derivatives of the propagator in Eq. ® one 
can use the following property 




£=0 


L 

/[+...+///—/ 


/ 



(Bl) 

Here the sum extends over all A-tuples of non¬ 

negative integers with h = I- 
Eor the derivatives we substitute 


de' 


u^{e{l+e)) 


E=0 


e'u^ 



(B2) 
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Appendix C: Numerical procedure 

Eirst, we construct a generic composite sequence of the 
form as shown in Eqs. © or ©. We calculate the derivatives 
from Eqs. ® or (fThl) using the identities (IBII) and (IB2I) . Then 
we proceed with a numerical minimization of the quantity 


Appendix A: Useful identities 

Erom the identity obtain 

T['(0/2)t/(0)Tj-(-0/2) = t/^(0), (Al) 

where E,(0) = and i denotes a certain qubit. Eor i = 2 

we have 


Uij, (0) = = cos 01 + isin OOxOij,. (A2) 


n 


^=E 

/=! 


a' 

de^ 


t/W(©)-t/(0) 


£=0 


(Cl) 


where we use Newton’s gradient-based method to determine 
the variables 0*. yielding D = 0. To minimize the number of 
the CPHASE gates t/^j, we start from a small number 
N, which is gradually increased, until we reach a solution to 
D = 0. Because we use a local optimization algorithm, we it¬ 
eratively pick the initial values of the variables using a Monte- 
Carlo scheme. 
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